Luttinger liquid physics from infinite-system DMRG 
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We study one-dimensional spinless fermions at zero and finite temperature T using the density 
matrix renormalization group. We consider nearest as well as next-nearest neighbor interactions; 
the latter render the system inaccessible by a Bethe ansatz treatment. Using an infinite-system 
alogrithm we demonstrate the emergence of Luttinger liquid physics at low energies for a variety of 
static correlation functions as well as for thermodynamic properties. The characteristic power law 
suppression of the momentum distribution n(k) function at T — can be directly observed over 
several orders of magnitude. At finite temperature, we show that n(k) obeys a scaling relation. The 
Luttinger liquid parameter and the renormalized Fermi velocity can be extracted from the density 
response function, the specific heat, and/or the susceptibility without the need to carry out any 
finite-size analysis. We illustrate that the energy scale below which Luttinger liquid power laws 
manifest vanishes as the half-filled system is driven into a gapped phase by large interactions. 

PACS numbers: 71.10.Fd, 71.10.Pm, 71.27.+a 



I. INTRODUCTION 

The concept of a Luttinger liquid (LL) is one of the 
most striking manifestations of strong electronic inter- 
actions in condensed matter theory^— It asserts that 
the low-energy excitations of a large class of gapless one- 
dimensional systems are not fermionic quasiparticles but 
collective (bosonic) modes; as a consequence, correla- 
tion functions exhibit anomalous power laws in space and 
time with interaction-dependent exponents. It is one of 
the hallmarks of a Luttinger liquid that the fermionic 
momentum distribution function n(k) exhibits a zcro- 
tcmperature power law singularity near the Fermi mo- 
mentum instead of a Fermi liquid discontinuity. 

The Luttinger liquid phenomenology is usually estab- 
lished by a two-step protocol. One first shows that un- 
der some rather restrictive assumptions - most impor- 
tantly of a linear dispersion ±v F \k — kp\ near the Fermi 
momentum hp and about the fc-dcpcndcnce of the two- 
particle potential g2,4,(k) - the low-energy limit of in- 
teracting fermions in one dimension is described by the 
Tomonaga-Luttinger (TL) model. In absence of the spin 
degree of freedom (on which we focus in the work), the 
TL Hamiltonian reads 
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where b n are free bosons, k n = 2?m/L, and L denotes the 
system size. The low-energy equilibrium physics of iJxL 
is determined by two scalars, the LL parameter K as well 
as the renormalized Fermi velocity vp, which depend on 
92.i(k = 0); introducing canonical fields LT and O yields 
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i?TL can be easily diagonalized, and after expressing 
fermionic annihilation operators in terms of the bosons 
b n one can derive analytical expressions for fermionic cor- 
relations functions at low energies. In a second step one 
needs to investigate about the validity of the initial as- 
sumptions that allowed to mathematically map interact- 
ing fermions to -Jztl- This is usually done by adding 
terms to the TL Hamiltonian which (partially) account 
for band curvature or physical choices of 52,4 (k). One 
can then employ renormalization group (RG) arguments 
to show that for a very general class of models those 
terms are RG irrelevant and that the low-energy physics 
is indeed described by a TL fixed point theory with renor- 
malized K and vp£~— 

Once a given Hamiltonian (one that might actually be 
motivated by an experimental setup and thus features 
a nonlinear dispersion, a finite bandwidth, and a phys- 
ical two-particle interaction) was shown to fall into the 
LL universality class, the behavior of thermodynamic ob- 
servables or correlation functions at low energies (or long 
length scales) is in principle known; it is universally de- 
termined by the expressions derived within the TL model 
if the corresponding values for K and vp are plugged inj 7 - 
However, there are only few quantitative verifications of 
Luttinger liquid physics directly for a microscopic the- 
ory, i.e. without resorting to the RG analysis within the 
TL fixed point model^ This particularly holds on the 
level of correlation functions - the power law singularity 
of the momentum distribution n{k), which is one of the 
key hallmarks of the LL phenomenology, has not been 
observed convincinglyJ^r— But even after a microscopic 
model was unambigously shown to exhibit LL physics 
at low energies, other important questions remain: (i) 
What is the energy (or temperature) scale below which 
LL behavior manifests? What is the physics away from 
the low-energy (zero-temperature) limit? (ii) How can K 
and Vp, which are non- universal functions of the model 
parameters, be obtained? (iii) How can the phase dia- 
gram be established? What happens to LL power laws if 
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the system is driven out of the LL phase? 

It is the first goal of this paper to address the above 
issues - most importantly to unambigously demonstrate 
the power law singularity of n(k) - for a lattice model 
of spinless fermions featuring nearest- as well as next- 
nearest neighbor interactions: 
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The ground state phase diagram for A2 = can be 
computed analytically via Bethe ansatzi 16 ' 17 The half- 
filled system is a Luttinger liquid if |A| < 1. A gap 
opens for large interactions as charge density wave or- 
der (A > 1) or phase separation (A < —1) develops. 
Away from half filling, the system is a LL for any re- 
pulsive A > 0. In presence of next-nearest neighbor 
interactions, no exact solution is known; it was shown 
numerically that the LL phase remains stable for small 
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One generally needs to resort to numerics (or ap- 



proximate analytics) if one wants to calculate correlation 
functions which are hard to obtain by Bethe ansatzi 18 i 19 
The density matrix renormalization group^ 3 . (DMRG) for 
finite systems was employed to investigate the fermionic 
momentum distribution^ 1 ^— the equal-time density re- 
sponse functio n 12 ! 24 " — as well as the density of states 
at zero temperature^2r— The results are consistent with 
the prediction from Luttinger liquid theory. The inverse 
system size, however, serves as an infrared cutoff, and 
despite substantial effort - systems of up to L = 3200 
sites were studied in Ref. _ very low energy (or long 
length) scales could not be accessed. To the best of our 
knowledge, fundamental LL power laws have not been ob- 
served unambigously over several orders of magnitude for 
any microscopic model within a numerically-exact frame- 
work. 

It is the second goal of this paper to propose the 
use of an infinite-system DMRG algorith m 23 i 31 i 32 (and 
to demonstrate its power) if one aims at studying Lut- 
tinger liquid physics at zero or finite temperature. Finite- 
size effects are absent by construction, and numerics are 
essentially controlled by one parameter \ (the dimen- 
sion of the DMRG block Hilbert spaces) only; as a side 
effect, a calculation for an infinite system is computa- 
tionally cheaper. We compute the single-particle Green 
function G(x) and its Fourier transform, the momen- 
tum distribution function n(k). At T = 0, the anoma- 
lous decay of G(x) ~ l/x 1+a , or equivalently the power 
law singularity of n(k) near the Fermi momentum Hf, 
\n(k) — 0.5 1 ~ \k — kp\ a , can be observed over several 
orders of magnitude; for A2 = 0, the exponent a is 
in perfect agreement with the TL formula 1 -! 2 - if K from 
Bethe ansat a 16 ' 17 is plugged in. At finite temperatures, 
we show that n(k) fulfills a simple scaling relation. We 



demonstrate that the energy scale on which the power 
law in n{k) manifests vanishes when large A or A2 drive 
the half-filled system into a gapped phase; the behavior 
of G{x) and n(k) in the latter is discussed. We more- 
over compute the equal-time density response function 
Cnn(x or k) at T — as well as the specific heat cy(T) 
and susceptibility £(T); from these quantities, the LL pa- 
rameters K and vf can be obtained without the need to 
perform finite-size or low-temperature extrapolations. 



II. INFINITE-SYSTEM DMRG 

We study the equilibrium properties of Eq. ([3]) at zero 
and finite temperature using a density matrix renormal- 
ization group algorith m 33 ' 34 for infinite systemsi 23 ' Any 
translationally-invariant state can be expressed in terms 
of a matrix product stated 



(4) 

where \a) denotes a local product basis, characterized e.g. 
by occupation numbers or a spin-1/2 degree of freedom 
after applying a Jordan- Wigner transformation. The ma- 
trices A^ 1 --- N ^ rTj are associated with a unit cell of size 
./V (see below). At zero temperature, we determine the 
ground state |^o) by applying an imaginary time evolu- 
tion exp(— tH) to a random initial state until the energy 
has converged to ten relevant digits. For T > 0, one first 
needs to purify the thermal density matrix px by intro- 
ducing an auxiliary Hilbert space Q: pr = Tr q\^t)(^t\- 
This is analytically possible only at T = 00 where 



I Pj, 00 3 Pj-oo 



(5) 



However, Y$>t) can be obtained from \%>oo) by applying 
an imaginary time evolution, |^t) = e~ g /( 2T ) l^oo) ' 36 ' 37 
Any static correlation function (on which we focus in this 
work) at zero or finite temperature can thus be exactly 
expressed as 



G AB (j u j 2 ) = (A h B h ) T = 



{* T \A n B j2 \^ T ) 



(6) 



For an orthogonal matrix product state,— calculating the 
thermodynamic limit expectation value (^x\Aj 1 Bj 2 \^>t) 
involves only a finite number of sites <jj S [ji, j'2]. 

In order to implement the above procedure, we em- 
ploy a standard time-dependent DMRG algorithm.— ~— 
After factorizing the evolution operators cxp(— ctH) us- 
ing a second order Trotter decomposition, they can be 
successively applied to Eq. (j4}. At each step, singular- 
value decompositions are carried out to update (some of 
the) matrices A^-" 1 ^"* within the unit cell. Since strict 
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FIG. 1. (Color online) Zero-temperature single-particle Green function G(x) = {Cj +x Cj)o and its Fourier transform, the 
fermionic momentum distribution function n(k), for a half-filled [quarter-filled in (b,iii)] lattice of spinless fermions exhibiting 
a nearest neighbor interaction A. Solid lines show DMRG data, dashed lines the asymptotic power laws G(x) ~ l/x 1+a( - K ^ 
and \n(k) — 0.5 1 ~ \k — k F \ expected for a Luttinger liquid. The LL parameter K is taken from Bethe ansatz i 16 ' 17 We 
use an infinite-system DMRG algorithm where the 'block Hilbert space dimension' \ i s the only numerical controll parameter. 
Finite-size effects are absent, and the LL power laws can be observed over several orders of magnitude. As the system is driven 
towards a gapped phase for A — > 1, the momentum scale on which they manifest vanishes. 



translational symmetry is broken by the Trotter decom- 
position, the latter is of size N = 2 for A2 = and 
N = 6 otherwise (at T = 0; for T > 0, one has N = 4 or 
N = 12, respectively). At zero temperature, the matrix 
dimension \ is fixed and constitutes the numerical con- 
trol parameter. For T > 0, x is dynamically increased 
during the imaginary time evolution of \^t) such that 
at each step the sum of all squared discarded singular 
values is kept below a threshold value e; the DMRG ap- 
proximation at finite temperature thus becomes exact in 
the limit e — > 0. We emphasize that finite-size effects arc 
by contruction absent in our approach. 

III. MOMENTUM DISTRIBUTION FUNCTION 

The single-particle Green function and its Fourier 
transform (the fermionic momentum distribution func- 
tion) are defined as 

G(x) = (c] +xCj ) T , n(k) =J2e- ikx G(x) . (7) 

X 

Within the exactly-solvable Tomonaga-Luttinger model 
of Eq. ([1]), G(x) and n{k) feature the following asymp- 
totic zero-temperature power lawsi — 

G(x > 1) - , \n(k « k F ) - 0.5| ~ \k - k F \ a , 

(8) 

for parameters where a < 1 on which we focus in this 
work [for a > 1 a linear term dominates in n(fc); singu- 
larities appear in derivatives] . The anomalous dimension 
a is a function of the LL parameter K only: 



Exploiting universality, one expects that Eq. ((HJ) also gov- 
erns the long-distance (or low-energy) behavior of our 
spinless lattice fermions. However, G(x) is hard to obtain 
by Bethe ansatz, and previous DMRG approaches^— 
employed finite systems of L = 0(100) sites where po- 
tential power laws were cut off by the inverse system size 
1/L. To the best of our knowledge, nobody has suc- 
ceeded in quantitatively demonstrating Eq. © - which 
is a hallmark of LL physics - for a microscopic model 
using a numerically-exact method. 

Within our infinite- system algorithm, G(x) is directly 
accesible via Eq. ([5]), and n(k) can be calculated by car- 
rying out the Fourier sum until it has converged [on the 
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FIG. 2. (Color online) The same as in Figure [T] but in pres- 
ence of next-nearest neighbor interactions A2 which render 
the system nonintegrable. Dashed lines show the power law 
~ l/x 1+a ( K ^ with the LL parameter K determined from the 
density response function. 
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x (k-k F )/n 

FIG. 3. (Color online) The same as in Figure [1] but for large interactions that drive the half-filled system into a gapped phase. 
G(x) decays exponentially, and n(k) is linear for close to fc_F. X i s chosen such that the results are converged on the scale of 
this plot. 



scale of the corresponding plot; typically, 0(10000) terms 
are necessary] . This allows us to verify the power laws of 
Eq. (J8| and to investigate on which energy or tempera- 
ture scale they manifest. We establish a scaling relation 
for n(/c, T) and provide strong indications that the latter 
is universal in the Luttinger liquid sense. 

A. Zero temperature 

We show the zero-temperature single-particle Green 
function G(x) as well as its Fourier transform n(k) in Fig- 
ures!]] [21 and[3l They feature asymptotic power laws over 
several orders of magnitude for any interaction strength 
A, A2 and filling where the system is a Luttinger liquid 
[the filling is always one half except for Figure [Ub,iii) 
where it is one quarter]. For A2 = 0, the corresponding 
exponent is in perfect agreement with Eq. @ if K from 
Bethe ansatz, which for half filling is given by 



is plugged in (dashed lines in Figure [T]). More quantita- 
tively, a power law fit of n(k) at A = 0.8 [at A = 0.25] 
yields an anomalous dimension a = 0.109 [a = 0.011] 
as compared to the exact result a(K Bcthc ) = 0.1096 
[a(-ftTBcthc) = 0.011]. In presence of next-nearest neigh- 
bor interactions, the model is nonintcgrable and K is not 
known exactly. In Section IIVI we will demonstrate that 
it can be extracted accurately from the density response 
function at small momenta; if we determine K for the pa- 
rameters of Figure [2] the corresponding exponent a{K) 
is consistent (dashed lines in Figure [2]). 

For small interactions, the power law decay of n(k) 
is cut off in the ultraviolet only by the bandwidth [see 
Figure [TJb,ii)]. At half filling, large interactions drive 
the half-filled system into a gapped phase: 4kp charge 
density wave order develops for A > 1, A2 = 0; if A2 > 
is increased at any fixed A, the system first enters a 



bond-ordered and eventually a 2kp charge density wave 
phased As one approaches the phase boundary, the scale 
on which the power law in n(k) manifests vanishes. This 
is illustrated for A2 = in Figure [T^b) ; note that even 
for A = 0.95 the power law is still obervable over two 
orders of magnitude in our infinite-system framework. 
In real space, the length scale above which G(x) features 
a power law becomes larger for A close to A = 1. The 
overall behavior is similar if A2 is increased at fixed A. 
Within the gapped phase, G{x) decays exponentially at 
large x (see Figure [3]). No LL power laws are visible at 
any intermediate scale. 

In our infinite-system DMRG approach at zero temper- 
ature, the matrices in Eq. Q which describe the ground 
state directly in the thermodynamic limit arc determined 
through an imaginary time evolution of a random ini- 
tial state . 23 ' 31 The key approximation results from the 
finite fixed matrix dimension \ in the state constructed 
by the algorithm (which leads to an effective correlation 
length that diverges as a power law in x)' 42 ' 43 However, 
no finite- x extrapolation is necessary for the results of 
this paper as reasonable matrix sizes already yield Lut- 
tinger liquid behavior over a broad range: Even moderate 
values x = 100 — 400 are sufficient to observe the power 
laws in G(x) or n{k) over several orders of magnitude 
(see Figured]). It should be emphasized again that x is 
the only numerical control parameter and that finite size 
effects are by construction absent in our approach. From 
a physical perspective, the use of an infinite-system al- 
gorithm seems reasonable from the beginning since in a 
finite system one would expect power laws to be cut off 
in the infrared by the inverse system size 1/L. Indeed, if 
we pragmatically restrict the Fourier sums in Eq. Q to 
O(100) sites - which is a typical value for a finite-system 
DMRG calculation 1 ^— - no power law behavior is vis- 
ible; it is moreover instructive to directly compare our 
Figure |Tfb) with Figure 3 of Ref . [l4| which was obtained 
using L = 66. Since it seems challenging to choose 1/L 
significantly smaller, the use of an infinite-system algo- 



5 









a=o.8, a 2 =o (a): 


10" 1 






■ 


10" 2 










; io 2 


t=i/ioo j 


10 -3 


: io 3 


V I 






' 10- 4 


~exp(-axf v ^ : 


\ y^-i/ioci^ 


10- 4 




50 100 150 


\T=1/20 \ 


10° 


10 1 


10 2 



O 
I 

ST 
c 








: A=0.iPX 

: scaling \ 
" relation 


\ — TL model 
: — . T=1/120 
- — T=1/140 
r — T=1/160 





0.1 1 
X 



0.001 



0.01 

(k-k F )/?i 



0.1 



FIG. 4. (Color online) Single-particle Green function and momentum distribution at finite temperature T. The LL power 
laws are cut off at long distances (small momenta) and replaced by an exponential (linear) decay. The infinite-system DMRG 
algorithm is controlled by the discarded weight e = 10" 11 during the imaginary time evolution of the density matrix from 
T = oo down to the physical temperature. In (b) , we compare our data (thin solid line) with the exact result (thin dashed line) 
for A = and T = 1/100 to demonstrate its accuracy. Left Inset to (b): Derivative dn(k)/dk\kp as a function of temperature. 
The solid line shows the expected TL power law ~ T"" 1 . Right Inset to (b): Rescaled variables y = T 1 ^ a dn(k)/d(Trk)\k F , 
x = vfQc — Isf) / {irT). Curves at different T collapse and are consistent with a scaling relation obtained for the TL model. 



rithm is advantageous if one aims at describing Luttingcr 
liquid physics of a microscopic model. 

The DMRG is frequently employed to compute the 
phase diagram of models that cannot be solved analyti- 
cally - e.g., the model at hand for A 2 > 0j2°-~— The type 
of order can be determined by considering various (e.g., 
charge density or bond) correlation functions. Careful 
finite-size analysis has to be carried out to accurately de- 
tect the phase boundary. This is not necessary within our 
approach. E.g., the transition out of the Luttinger liquid 
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FIG. 5. (Color online) Zero-temperature spin correlation 
function S(j) — (S^Sq) of a spin-1/2 Heisenberg chain with 
z-anisotropy A. Solid lines shows DMRG data, dashed lines 
a field theory result of Ref. [45l . The latter yields both the 
exponents and coefficients of an asymptotic expansion at 
large distances. Inset: S 2 (j) = [8S{j) - 8S(J + l)]/2 with 
5S(j) = S(j) - Si(j) and S&(j) = A/f being the leading 
term within field theory. The difference Sa(J) is governed by 
the leading alternating term. 



phase at A = 1, A 2 = can be easily identified from the 
behavior of G(x) using moderate values of x ~ 500 [see, 
e.g., the curve at A = 1.04 in Figure [3{a) which clearly 
shows that the system has entered a non-LL phase]. 



B. Finite temperature 

We now turn to finite temperatures. This is interesting 
since (a) LL power laws will be cut off in the infrared by 
T; thus, one needs to study on which temperature scale 
they are still observable for a given microscopic model 
(e.g. one that could be motivated by an experimental 
setup), (b) little is known about n(k, T) if both \k — kp\ ^ 
and T^O; due to the absence of any low energy scale 
one might expect a scaling relation for n(k,T), and (c) 
adressing those questions establishes the infinite-system 
DMRG algorithm as a convenient tool to study systems 
at nonzero temperature directly in the thermodynamic 
limit. 

Let us briefly recall some computational details: The 
single-particle Green function G{x) (or any other corre- 
lation function) can be calculated staightforwardly via 
Eq. (|6|), and the state \^t) which purifies the density 
matrix is obtained from the one at T = oo (where it 
is known) through an imaginary time evolution down to 
the physical temperature Ti 36 i 37 As mentioned, we again 
employ an infinite-system algorithm, i.e. work with a ma- 
trix product state that describes a state directly in the 
thermodynamic limiti 23 ' 31 Finite-size effects are thus by 
construction absent, and the computational effort is mod- 
erate since only N = 4 (for A 2 = 0; N = 12 for A 2 > 0) 
different matrices appear in Eq. @. Numerics are con- 
trolled by the discarded weight e during each Trotter step, 
which we typically fix to e = 10~ n (for a step size of 
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FIG. 6. (Color online) DMRG calculation of the zero-temperature instantaneous density response function Cnn(x) = 
{ c \+ x c 3+x c \ c j)o- Its Fourier transform Cjvjv(fc) exhibits a power law singularity Cjvjv(A; ~ 2k f) ~ (2kp — k) 2K ~ 1 for mo- 
menta near 2&f. For small k, Cjvjv(fc) is linear with a slope determined by the Luttinger liquid parameter K. The latter can 
thus be determined accurately through linear fits [see the comparison with Bethe ansatz results in (b)]. 



At = 0.01). A comparison with the analytic result at 
A = A 2 = and T = 1/100 shows that this is sufficient 
to obtain numerically exact data [see Figure HJb)]. 

Our results are shown in Figure [4] Power laws are 
cut off by the temperature; G(x) decays exponentially at 
large lengths, and n(k) is linear for momenta close to &f- 
For T < 1/100 and A = 0.8, the LL power law in n(k) 
is observable over about one order of magnitude in \k — 
kp\. Due to the absence of any low energy scale except 
for vp\k — kp\ and T, one also expects power laws in 
temperature. Indeed, for the TL model the fc-derivative 
satisfies^ 



dn(k) 



,1k 



(11) 



which is consistent with our data for the lattice model 
[see the left inset to Fig. [4fb)]. For \k — kp\ ^ and 
T ^ one can show that within the TL model n(k,T) 
has a simple scaling form 



dn(k) 
dk 



rpa-lp 



Vp(k — kp) 



ttT 



(12) 



Eq. (|T2"j) can be established (and the scaling function F 
can be computed) by expressing n(k,T) in terms of the 
local density of states which within the TL model was 
shown to have scaling form4i In the right inset to Figure 
2Jb) we illustrate that our DMRG data for the lattice 
fermions at small T and \k — kp\ agrees with Eq. (fT2")) and 
with the form of F (if in the latter both axes are scaled 
arbitrarily). This strongly indicates that the low-energy 
behavior of n{k, T) at finite temperature is universal for 
any model that falls into the Luttinger liquid universality 
class. 



C. Spin chain 

For the integrable model, we have shown that G(x) 
asymptotically decays as ~ l/x 1+a ^ K > if the LL param- 
eter K from Bethe ansatz is plugged in. We will now 
investigate if one can make a reasonable statement about 
subleading terms. To this end, we consider chain of spin- 
1/2 degrees of freedom S x ' y,z : 

H s = {S*S* +1 + Spy +1 + &S*S* +1 ) , (13) 



which can be mapped to Eq. ([3]) through a Jordan- 
Wigner transformation 



Si 



(14) 



For the longitudinal correlation function S(j) = (SjSq)o 
- which is the spin equivalent of the single particle Green 
function of Eq. ([7]) - an asymptotic expansion 



S(j) 



J 



A B , .. C 

— + (-iy— 



in a 



(15) 



was obtained by field theory,— and the exponents as well 
as the prefactors were computed. One might speculate 
that the string operators in Eq. (|14[) hinder a similar 
approach in the fermionic case. 

In Figure[5]we show that zero-temperature DMRG cal- 
culations for S(j) at A = 0.8 compare nicely with the 
field theory result^ The exponent of the leading term 
r\A is much smaller than rjB,c and thus dominates at large 
lengths. If we substract this term, SS(j) = S(j) — A/j VA , 
and compute the difference S2C7) = [$S(j)—5S(j + l)}/2, 
the leading alternating term (— l) J C/r]c can clearly be 
observed. The subleading homogeneous term B/j VB is 
consistent with the behavior of S^ij) = [SS(j) + SS(j + 
l)]/2 but cannot be detected with similar accuracy (be- 
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FIG. 7. (Color online) Thermodynamic quantities from 
infinite-system DMRG. The main panel shows the suscep- 
tibility £(T) = T- 1 J2 j (SfSo)t, = c] C j - 1/2; the inset 
displays the specific heat ey(T). The low-temperature be- 
havior is determined by the LL parameter K and renormal- 
ized Fermi velocity vf. At A 2 = 0, the dashed lines show 
the asymptotic behavior predicted by Bethe ansatz.— In the 
nonintegrable case, our approach provides a simple way to 
extract K and vf from £ and cy without the need to carry 
out finite-size (low-temperature) extrapolations. 



cause any slight derivation from the substracted field the- 
ory exponent i]a significantly modifies <S^~ but not the 
difference S2). 



IV. EXTRACTION OF LL PARAMETERS 

In this Section we discuss further static correlation 
functions as well as thermodynamic quantities. Our main 
focus will be to demonstrate how the Luttinger liquid pa- 
rameter K as well as the renormalized Fermi velocity vf 
- which fully characterize the TL model and thus LL 
physics - can be extracted without the need to carry out 
any finite-size or low-temperature extrapolations ! 25 ' 26 
For A2 = we achieve excellent agreement with the 
exact Bethe ansatz result. Infinite-system DMRG thus 
provides a reliable way to determine K and vp for a non- 
integrable model. 



A. Density response 

The two-particle analogue of G(x) is the instantaneous 
density response function 



C NN {x) 
C NN (k) = 



C j+x C j+x C 



-ikx 



nn(x) 



(16) 



Within the Tomonaga-Luttingcr model, it exhibits the 
following characteristic zero-temperature behavior,i2 



C N n {x) 



K 



2ir 2 x 2 



C 



cos(2Icfx) 



r 2K 



(17) 



where C is a nonuniversal constant. In momentum space 
this translates to 



dC NN {k) 
dk 



K/{2%) 
~ (2k F - 



k) 



2K-2 





2k F 



(18) 



12.24-26 



at 



DMRG data for the spinless lattice fermionS' 
half filling is shown in Figure |H1 We reproduce the 
power law singularity for momenta close to 2kp; a power 
law fit at A = 0.8, A2 = yields an exponent of 
—0.77 which is in good agreement with the exact result 
2i\Bcthc - 2« —0.74. The slope at small k is determined 
by the LL parameter K which can thus be determined 
accurately from simple linear fits. Even for moderate 
values of x — 200 (which is the only numerical con- 
trol parameter) , K agrees perfectly with the exact Bethe 
ansatz result Eq. (fTOj) for A2 = [see the comparisons at 
A = 0.25, 0.5, 0.8 in FigureEJb)]. Infinite-system DMRG 
thus provides a reliable way to compute K for noninte- 
grable models. 

The LL parameter K was previously calculated from 
the density response at small momenta using DMRG for 
finite systemsi^ 5 . It turned out that using L k, 100 — 400 
sites is sufficient to extrapolate to the thermodynamic 
limit, which seems reasonable since the first term in 
Eq. (fT7)) falls off rapidly. The advantage of our approach 
is that finite-size extrapolations are superfious and that 
the computational effort is minimal since only a single 
unit cell (and not L — 400 sites) needs to be considered. 



B. Thermodynamic quantities 

We finally consider two thermodynamic quantities: the 
specific heatii 2 - 



c v (T) = 



d(H) T t^o T 



dT 



Vf 



(19) 



as well as the susceptibility (we use the terminology 'sus- 
ceptibility' following previous works on the Hciscnberg 
chain)ii 



1/2 



1/2 



Vf 



(20) 

Both the LL parameter K and the renormalized Fermi 
velocity vf can be determined from the low-temperature 
behavior of cy/cy and £/f°, where the ratio with the 
noninteracting (known) values Cy and £° is taken so 
that nonuniversal prefactors drop out. In Figure [71 we 
show cy and £ in comparison with the Bethe ansatz 
asymptotes^ for A2 = 0. The low-energy regime can be 
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reached easily, and the agreement with the exact result is 
quantitative. Thus, the infinite-system DMRG algorithm 
at nonzero temperature provides an accurate means to 
determine K and Vp for nonintegrable models without 
the need to carry out low-temperature extrapolations^ 

V. CONCLUSION AND OUTLOOK 

We studied zero- and finite-temperature Luttingcr liq- 
uid physics of spinless lattice fermions exhibiting nearest 
and (nonintegrable) next-nearest neighbor interactions 
terms A and A2. Employing an infinite-system density 
matrix renormalization group algorithm (and carrying 
out Fourier sums up to <~ 10000 sites), we unambigously 
showed that the momentum distribution function n{k) at 
T = exhibits the characteristic power law suppression 
near the Fermi momentum. For A 2 = 0, the associated 
exponent is in perfect agreement with the one predicted 
for the Tomonaga-Luttinger model if the corresponding 
LL parameter K is taken from Bcthe ansatz. The scale on 
which this power law is visible vanishes when large inter- 
actions drive the half-filled system into a gapped phase. 



Finite temperatures cut off power laws in the infrared, 
and n(k,T) obeys a scaling relation. We illustrated that 
the LL parameter K and renormalizcd Fermi velocity vf 
are easily extracted from the instantaneous density re- 
sponse at small momenta, from the specific heat, or from 
the susceptibility; the transition into a non-LL phase can 
be easily detected. Both is of special importance for mod- 
els which cannot be solved analytically. Any finite-size 
analysis is superflous, and our algorithm is computation- 
ally cheap since only a single unit cell needs to be con- 
sidered. This overall illustrates that employing infinite- 
system DMRG is advantageous if one aims at investi- 
gating Luttinger liquid physics of microscopic models at 
low energies or temperatures. As a next logical step it 
would be interesting to calculate time- dependent corre- 
lation functions, and their Fourier transform, the local 
density of states. 
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